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ABSTRACT 

We have discovered that the brown dwarf WISE J014656.66+423410.0 is a close bi¬ 
nary (0 / /0875±0'. / 0021, 0.93±g;i| AU) from Keck laser guide star adaptive optics imaging. 
Our photometry for this system reveals that both components are less luminous than 
those in any known substellar binary. Combining a new integrated-light spectrum (T9p) 
and resolved Y JH -band photometry from Keck allows us to perform spectral decompo¬ 
sition and assign component types of T9 and Y0. Many of the unusual features in the 
spectrum might be explained by high surface gravity: Y-band peak broadened to the 
blue; J-band peak broadened to the red; R-band peak shifted slightly to the red; and 
red Y — J colors. Interestingly, the very low component luminosities imply that the T9 
primary is unexpectedly cold (T e g = 345 ± 45 K assuming an age of 10 Gyr), making it 
«100K cooler than any other late-T dwarf and comparable to Y dwarfs. One intriguing 
explanation for this apparent discrepancy is that the J- and R-band spectral features 
that trigger the transition from T to Y spectral types are highly gravity-dependent. 
This can be tested directly in the very near future by orbit monitoring. We constrain 
the orbital period to be <10 yr by combining evolutionary model-based mass estimates 
for the components (~12-21 Mj up , lcr at 10 Gyr) with a statistical constraint on the 
semimajor axis (<1.3 AU). Such a period is shorter than any other known T/Y tran¬ 
sition binary, meaning that WISE J0146+4234AB will likely yield a dynamical mass 
within the next few years. 
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1. Introduction 

The last few years have seen a rapid increase in the discovery of brown dwarfs significantly 
cooler than any previously known objects in the solar neighborhood (Luhman et al. 2011; Liu et al. 
2011; Kirkpatrick et al. 2011, 2012; Luhman 2014). Cushing et al. (2011) presented the discovery of 
several brown dwarfs from Wide-field Infrared Survey Explorer (WISE) data, and they distinguished 
these objects from warmer brown dwarfs primarily using spectral features in the J and H bands 
and proposed a new “Y” spectral class for the sample. Dupuy & Kraus (2013) presented the first 
comprehensive sample of parallaxes for the Y dwarfs and derived bolometric luminosities showing 
that even objects with quite different near-infrared spectra (and thus very different T e ff estimates 
from model atmospheres) in fact had indistinguishable temperatures of ~400-450 K. Assuming ages 
of 1-5 Gyr, typical for field brown dwarfs, evolutionary models predict masses of 6-20 Mj up for the 
Y dwarfs. Thus, Y dwarfs could be free-floating planetary-mass (<13Mj up ) objects, but this can 
only be tested directly with dynamical mass measurements from binaries. However, none of the 
binaries yet discovered at the T/Y transition (Gelino et al. 2011; Liu et al. 2011, 2012) have short 
enough estimated periods (<20yr) to enable orbit determinations in the near future. 

We present here the discovery of a tight (0709) binary at the Y dwarf boundary from our 
ongoing Keck laser guide star adaptive optics (LGS AO) imaging survey. Mace et al. (2013a) 
discovered the unresolved source WISE J014656.66+423410.0, hereinafter WISE J0146+4234, and 
assigned a spectral type of Y0 based on Keck/NIRSPEC spectroscopy in J and H bands (also see 
Kirkpatrick et al. 2012). Using a new, higher signal-to-noise ratio (S/N) Gemini/GNIRS spectrum 
covering 0.9-2.4 pm we derive an integrated-light spectral type of T9p for WISE J0146+4234. 
Adopting the parallax distance of 10.6^'g pc from Beichman et al. (2014), the projected separation 
of WISE J0146+4234AB is ~0.9 AU, which would make it tighter than any previously known L or 
T dwarf visual binary at discovery. 


2. Observations 
2.1. Keck/NIRC2 LGS AO 

We first observed WISE J0146+4234 during a period of exceptionally good seeing on 2012 Sep 7 UT 
with the facility near-infrared camera NIRC2 with the LGS AO system at the Keck II telescope 
(Wizinowich et al. 2006; van Dam et al. 2006). According to the differential image motion monitor 
(DIMM) at the Canada-France-Hawaii Telescope, the seeing in the visible was 0 / . / 35 on Mauna Kea 
during our observations. Compared to all archived seeing measurements during 2012, such favor¬ 
able conditions occur <5% of the time that the DIMM is active. During our dither sequences we 
kept the LGS centered in NIRC2’s narrow camera field of view. The wavefront sensor recorded 
an LGS flux equivalent to a L ~ 9.8-10.2 mag star and a lower bandwidth sensor monitored the 
star 2MASS J01465144+4234388 (R = 15.3 mag, 65" away from the target) in order to perform 
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tip-tilt corrections. We later observed WISE J0146+4234 in the same setup on 2012 Oct 8 and 
2013 Oct 22 UT, during times of unusually good seeing (« 0 / . , 30~0 / . / 40 and 0 / .'45-0 / /55, respectively). 

We reduced our LGS AO data in the same fashion as our previous work (Liu et al. 2008; Dupuy 
et al. 2009b,c, 2010). Reduced images from each data set are shown in Figure 1. We measured 
binary parameters by fitting three-component Gaussians to individual images and determined the 
uncertainties in these parameters from the rms scatter among each data set. To correct for non¬ 
linear distortions in NIRC2 we used the calibration of Yelda et al. (2010), from which we also adopt 
linear terms of 9.952 ± 0.002 mas pixel” 1 for the pixel scale and +0?252 ± 0?009 for the detector’s 
+y -axis orientation relative to North. Table 1 presents the binary separation, position angle (P.A.), 
and flux ratio that we derived from each data set, along with the airmass and AO-corrected FWHM 
of our images. We note that the astrometry measured in different bandpasses at a given epoch agree 
well within the errors and have reasonable x 2 ) supporting our adopted astrometric uncertainties. 


2.2. Gemini/NIRI Photometry 

WISE J0146+4234 had no previously published photometry on the Mauna Kea Observatories 
(MKO) system, and the 2MASS system photometry from Kirkpatrick et al. (2012) has very large 
errors, so we obtained new MKO photometry using the facility camera NIRI at the Gemini-North 
Telescope (Hodapp et al. 2003) on 2013 Jan 10 and 2013 Jan 12 UT. All data were taken on photo¬ 
metric nights, with seeing typically around (//8 and UKIRT FS 6 used for photometric calibration. 
Our observations at Y, J, H , CH^s, and K bands are summarized in Table 2. We note that our 
photometry is highly inconsistent with the values presented by Kirkpatrick et al. (2012), even after 
accounting for the conversion from 2MASS to MKO systems. Jmko — J 2 MASS = —0.260 mag and 
Hmko — JJ 2 MASS = 0.074 mag according to synthesized photometry from our GNIRS spectrum 
(Section 2.3). Our photometry is 1.6mag (6. Oct) fainter at Jmko band and 2.5mag (9.3cr) fainter 
at JJmko band. The cause of this discrepancy is not clear, but we note that Leggett et al. (2013) 
reported similar cases for three Y dwarfs that they found were significantly fainter (0.5-1.0 mag) 
than published values from Cushing et al. (2011) and Kirkpatrick et al. (2012). We suspect large 
systematic errors in the photometry from Kirkpatrick et al. (2012) for WISE J0146+4234 because 
colors based on those values were noted by Lodieu et al. (2013) to be extremely unusual (z — J 
and z — H colors >2mag redder than objects of similar spectral type), whereas our J and H 
magnitudes resolve this discrepancy. Beichman et al. (2014) have also reported iL-band photom¬ 
etry (Hmko = 20.91 ± 0.21 mag) that agrees with our higher precision Gemini measurement. We 
therefore use our photometry in the following analysis. 
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2.3. Gemini/GNIRS Spectroscopy 

We obtained a full 0.9-2.4 /xm spectrum of WISE J0146+4234 using the GNIRS spectrograph 
(Elias et al. 2006) because the published spectrum from Mace et al. (2013a) had limited wavelength 
coverage and low S/N. We used the 32 lines mm -1 grating, 0 / . , 15 pixel -1 camera, and 0 / . , 675 slit to 
achieve a resolving power of R ~ 750 over a total integration time of 8400 s (Table 2). We reduced 
our data in a similar fashion to our previous work (e.g., Leggett et al. 2014a), and the flat fielded, sky 
subtracted, and rectified integrated-light spectrum of WISE J0146+4234AB is shown in Figure 2. 
We flux calibrated the spectrum using our NIRI photometry at Y, J, and H bands (Table 2). To 
account for modest discrepancies between our synthesized and measured YJH colors, we added 
0.2 mag errors in quadrature to our synthesized photometry to achieve p(x 2 ) = 0.5 when computing 
the optimal scaling factor and its error. 

From our calibrated spectrum we computed AT-band photometry as check against our very 
low S/N NIRI photometry (AT = 22.4 ± 0.4mag). After accounting for the uncertainty in the 
calibration, we synthesized K = 21.31 ± 0.24 mag. This is inconsistent with our NIRI photometry 
at 2.3a. Both our photometry and the A'-band portion of our GNIRS spectrum are very low S/N, 
and our synthesized photometry may be subject to systematic errors associated with extrapolating 
our flux calibration from Y JH to K. Therefore, in the following analysis we adopt an average value 
of K = 21.75T0.25 mag that is consistent at 1.3(7 with both measured and synthesized photometry. 


3. Results 
3.1. Companionship 

From our resolved photometry alone, companionship is highly likely. The blue near-infrared 
colors of the companion WISE J0146+4234B are very similar to WISE J0146+4234A, implying that 
it not likely to be a background giant. The difference in CH±s — H color of 0.29 ± 0.15 mag between 
the two components implies that they have similar levels of methane absorption and therefore 
are both late-type T or Y dwarfs. According to Liu et al. (2012), the probability of finding such 
an unassociated late-T dwarf in the entire l(/. , 2 x 10 / . , 2 field of view of NIRC2’s narrow camera is 
1.6 x 10 -6 , thus the likelihood that the two components of WISE J0146+4234AB are not physically 
bound is negligible. Moreover, the proper motion of (p a cos 6 , ps) = (—O'/441 ± 07013, —0 , /026 ± 
0 / /016) yr -1 and parallax of 0 / . , 094 ± O'/Old measured by Beichman et al. (2014) implies that an 
unassociated background star would move by 0 / . / 56 ± 0 / /02 between our first and last observations, 
which is strongly ruled out by our astrometry (Table 1). 
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3.2. Spectral Types 

To determine the integrated-light spectral type of WISE J0146+4234AB we directly compared 
our spectrum to standards from Cushing et al. (2011). For reference, we also computed commonly 
used spectral indices, which we report in Table 3. In determining spectral types we focus on the 
J- and 17-band regions of the spectrum, leaving the somewhat unusual Y-band region to our later 
discussion in Section 3.5.1. 

The relative heights of the flux peaks in the J and H bands better match the T9 standard 
UGPS J072227.51—054031.2 and T9.5 dwarf WISEPC J014807.25—720258.7 than the Y0 standard 
WISEPA J173835.53+273258.9 (Figure 2). The J-band region, normalized to its peak flux, best 
matches the T9 standard UGPS J0722—0540 because the Y0 standard, and even the T9.5 dwarf 
WISE J0148—7202, have notably narrower flux peaks. This is reflected by the fact that J-band 
spectral indices W j (Warren et al. 2007) and J-narrow (Mace et al. 2013a) best match the T9 dwarfs 
from Mace et al. (2013a). However, the peak-normalized 17-band region seems to be intermediate 
between T9 and Y0 standards, best matching the T9.5 dwarf WISE J0148—7202. In terms of 
indices, this manifests as an NH 3 — H index (Delorme et al. 2008) that best matches T9.5 dwarfs. 
We also note that the wavelength at the peak of the Jl-band flux is shifted slightly to the red 
by 0.006/rm, a feature that has only previously been observed in WISEPC J140518.40+553421.4, 
which was typed as Y0(pec?) by Cushing et al. (2011). No two of these three traits (T9-like J band, 
T9.5-like H band, and shifted 17-band peak) have been observed in any other object, making the 
integrated-light spectrum of WISE J0146+4234AB particularly unusual. The two methane indices 
(defined by Burgasser et al. 2006b) are broadly consistent with a late-T or Y spectral type. The 
CH 4 — J index best matches T9.5 dwarfs but is also consistent with Y0 dwarfs. The CH 4 — H 
index is problematic because the numerator has almost no flux in it by the late-T dwarfs. Thus, 
even though our spectrum has sufficient S/N to accurately measure this index, many comparison 
objects do not, so our value ends up being consistent within the rms of all T8.5-Y0 spectral type 
bins. 


Overall, we find evidence for both T9 and T9.5 spectral types by visually comparing to spectral 
standards and examining spectral indices. Most published work has focused on J band in deter¬ 
mining spectral types at the T/Y transition, and our J band spectrum leans more toward a type of 
T9. Therefore, we assign a type of T9p to the integrated-light spectrum of WISE J0146+4234AB. 
In this case, the peculiar designation refers to both the 0.006 fi m shift in the 17-band peak, as 
originally suggested by Cushing et al. (2011) for WISE J1405+5534, and the T9.5-like H -band 
spectrum. 

To determine spectral types for the individual components of WISE J0146+4234AB, we per¬ 
formed spectral decomposition following the method outlined in Section 5.2 of Dupuy & Liu (2012). 
Briefly, we consider all possible pairs from a library of template spectra and find the optimal scale 
factors needed for a pair of templates to best match the observed integrated-light spectrum. We 
then compute synthetic relative photometry for these pairings and compute the y 2 of these val- 
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ues compared to the measured flux ratios from our Keck LGS AO imaging. Unlike in Dupuy 
& Liu (2012), here we have a much smaller library of template spectra given the number of ob¬ 
jects later than T 8 with high quality near-IR spectra. Our library comprises four T8.5 dwarfs 
(ULAS J003402.77—005206.7, WISEPC J121756.91+162640.2A, ULAS J133553.45+113005.2, and 
Wolf 940B), two T9 dwarfs (UGPS J0722—0540 and CFBDSIR J145829+101343AB), one T9.5 
dwarf (WISE J0148-7202), five Y0 dwarfs (WISEPA J041022.71+150248.5, WISE J1217+1626B, 
WISE J1405+5534, WISE J1738+2732, and WISEPC J205628.90+145953.3), and one Y0.5 
(WISEPA J154151.66—225025.2). Spectra for these objects were published by Warren et al. (2007), 
Burningham et al. (2008, 2009), Cushing et al. (2011), and Leggett et al. (2014a). Another dif¬ 
ference in our analysis compared to Dupuy & Liu (2012) is the wavelength range tested. Because 
many published spectra do not extend into Y or K band, we use only the 1.15-1.85 pm region for 
our fitting. 

The best matching template pair, both in terms of fitting our integrated-light spectrum and 
flux ratios at J, H , and CH^s bands, is UGPS J0722—0540 (T9) and WISE J1405+5534 (YOp). 
A primary spectral type of T9 is consistent with the fact that this brighter component should 
dominate the integrated-light spectrum given our Keck flux ratios of ~1.0mag in J and H bands. 
It is not surprising that WISE J1405+5534 provides the best matching secondary component be¬ 
cause of some of the features it shares with our integrated-light spectrum (shifted fL-band peak 
and large CH 4 — H index). Figure 3 shows this best-fit template pairing. A number of other 
template pairings give both a good fit to the spectrum and reasonable flux ratios. These either 
use UGPS J0722—0540 (T9) or a T8.5 (ULAS J0034—0052 or WISE J1217+1626A) as the pri¬ 
mary and a Y0 (WISE J1217+1626B, WISE J1738+2732, or WISE J2056+1459) or the T9.5 
WISE J0148—7202 as the secondary. Therefore, we assign component types corresponding to our 
best-fit template pair (T9+Y0) with uncertainties of ±0.5 subtypes in each. Given some of the 
unusual features seen in integrated light, we note that one or both components may also be peculiar 
despite our not typing them as such. 


3.3. Near- and Mid-Infrared Colors and Magnitudes 

Our new photometry and spectral analysis allow us to compare the colors of WISE J0146±4234AB 
to other late-T and Y dwarfs (Figure 4). The Y JH colors of both components and the integrated- 
light J — K color are mostly typical of other T8.5-Y0 objects. The color in which the components 
appear most unusual is Y — J, where the Y0 secondary’s color is 0.80 ±0.20 mag. This would not be 
unusual for a late-T dwarf, but it is on the extreme red end of colors for Y dwarfs where only one 
comparably red object is known, the Y0 dwarf WISE J073444.02—715744.0 ( Y—J = 0.97±0.07 mag; 
Leggett et al. 2014b). Moreover, both components of WISE J0146±4234AB are distinctive in that 
they are the reddest known objects in Y — J at faint absolute magnitude (My > 21 mag). Their 
J — H colors place them on the blue edge of the T/Y sequence in color-magnitude diagrams, and 
their integrated-light J — K color of —1.06 ± 0.26 mag is normal for objects of comparable absolute 
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magnitude but would be relatively blue for a late-T dwarf. 

On mid-infrared color-magnitude diagrams, the integrated-light photometry of WISE J0146+4234AB 
follows the sequence of late-T and Y dwarfs, but its location is unusual given the integrated- 
light spectral type of T9. WISE J0146+4234AB’s integrated-light color [3.6] — [4.5] = 2.42 ± 

0.07 mag is the reddest of any known T9 dwarf, with the next reddest objects being the T8.5+Y0 
binary WISE J1217+1626AB (T9 in integrated-light with [3.6] — [4.5] = 2.33 ± 0.03 mag) fol¬ 
lowed by WISE ,1033515.01+431045.1 (T9) and WISEPA J213456.73—713743.6 (T9p) both with 
[3.6] — [4.5] = 2.22 ± 0.04 mag. The only late-T dwarf that is as red as WISE J0146+4234AB 
is WISE J081117.81—805141.3 (T9.5:, [3.6] - [4.5] = 2.42 ± 0.07mag), but like other T9 and 
T9.5 dwarfs it is ~1.0—1.5 mag brighter in [3.6]-band absolute magnitude than the combined 
light of WISE J0146+4234AB. The situation is similar but less extreme in J — W2 color ver¬ 
sus W2-band absolute magnitude. WISE J0146+4234AB is the reddest late-T dwarf (J — W2 = 

5.61 ± 0.10 mag) with the exceptions of WISEPC J232519.54—410534.9 (T9p; 5.64 + 0.04 mag) and 
WISE J2134—7137 (T9p; 5.86 ± 0.11 mag). Both of these were typed as peculiar by Kirkpatrick 
et al. (2011) due to excess flux at Y band and less flux at K band, similar to features seen in 
the integrated-light spectrum of WISE J0146+4234AB. However, unlike these T9p objects, the 
integrated-light VP2-band flux of WISE J0146+4234AB is 0.7-0.8 mag fainter, implying that the 
individual components are >lmag fainter. 

Both components of WISE J0146+4234AB have very faint absolute magnitudes given their 
spectral types. WISE J0146+4234A is the faintest T9 dwarf known in any near- or mid-infrared 
bandpass, and WISE J0146+4234B is as faint or fainter than WISE J1405+5534 (YOp) in near- 
infrared magnitudes and could be the faintest in the mid-infrared, too, depending on the binary 
flux ratio at those wavelengths. WISE J0146+4234A is the only late-T dwarf wholly overlapping 
in colors and magnitudes with the Y dwarfs. 


3.4. Bolometric Magnitudes 

Dupuy & Kraus (2013) showed that summing fluxes across the near- and mid-infrared can 
produce an estimate for the bolometric magnitudes (mt, 0 i) of late-T and Y dwarfs that is only weakly 
dependent on the assumed model atmospheres. We applied this “super-magnitude” approach to 
the integrated-light photometry of WISE J0146+4234AB and found mboi = 21.58 ± 0.12 mag. To 
apportion this between the two individual components we employed our resolved near-infrared 
photometry. For our purposes, only the relative near-infrared super-magnitude between the two 
components is relevant, and we find A myjH = 0.95 ± 0.05 mag. To convert this into a bolometric 
flux ratio requires an estimate of the bolometric correction for each component. Dupuy & Kraus 
(2013) report bolometric corrections of BC yjh = 1-6 ± 0.6 mag at T9 and 0.8 ± 0.6 mag at Y0, 
and we adopt these values for the primary and secondary, respectively. Thus, the bolometric flux 
ratio is calculated as Ambol = A myjH + BCyj^(YO) — BCyj#(T9) = 0.2 ± 0.8 mag, and we find 
bolometric magnitudes of mbol = 22.2 ±0.4 mag and 22.4 ±0.4 mag for the primary and secondary, 


respectively. 


We compare our apparent bolometric magnitude estimates to the typical absolute magnitudes 
of late-T and Y dwarfs in order to illustrate an unusual property of WISE J0146+4234AB. Using 
the Lboi values from Table S5 of Dupuy & Kraus (2013), we find that M^ 0 \ = 20.1 ± 0.5 mag 
for T9 and Mboi = 21.06 ± 0.25 mag for Y0, where the error bars represent the intrinsic rms 
scatter in bolometric magnitude among objects of the same spectral type. Assuming spectral types 
of T9+Y0 for WISE J0146+4234AB, the integrated-light absolute bolometric magnitude would be 
Afbol = 19.7±0.4 mag if they were typical of other known objects. However, this is 1.7 mag different 
from the bolometric magnitude of M^ 0 \ = 21.41 q3 mag derived from the photometry above and the 
parallactic distance of 10.6ljgpc from Beichman et al. (2014). Therefore WISE J0146±4234AB 
is significantly less luminous than other known objects of the same spectral type, which at least 
partly explains its unusual appearance on color-magnitude diagrams as described above. 

We consider the alternative that the parallax of 94 ± 14 mas from Beichman et al. (2014) is in 
error. They used multiple telescopes (WISE, Spitzer, and Keck) and bandpasses ( H , [3.6], [4.5]) 
to measure the parallax, which introduces the possibility of systematic errors in the astrometry. 
This is especially true in the case of a T9±Y0 binary where the flux ratio could vary considerably 
from the near- to mid-infrared causing unaccounted for shifts in center-of-light measurements at 
different epochs. However, for the components of WISE J0146±4234AB to have normal magnitudes 
for a T9±Y0 dwarf pair would require a distance of 26^g pc, where the error bars account for the 
uncertainty in the mboi and scatter in M^ 0 \ as a function of spectral type (Dupuy & Kraus 2013), 
and thus a parallax of «40 mas. This would require photocenter shifts of more than half the binary 
separation, and such shifts would have to be correlated with the parallax factor. Figure 17 of 
Beichman et al. (2014) shows that there is both Keck and Spitzer data on both sides of the ellipse, 
so a large error in the parallax due to mixing bandpasses seems unlikely. 

We also consider whether orbital motion could have impacted their parallax. Because proper 
motion accommodates for any linear orbital motion, only acceleration is a concern, and this usually 
shows up in the residuals of the fit. Beichman et al. (2014) report x 2 = 23.0 (27 degrees of 
freedom), so there are no significant residuals due to orbital motion, but such acceleration could 
in principle be aligned with the parallax motion by chance. To test this possibility, we fitted 
just the Beichman et al. (2014) astrometry spanning 2012 October to 2013 November, which is 
contemporaneous with our Keck data that show only a small amount of orbital motion, (d, 5) = 
(—4.3 ± 3.9, 3.9 ± 2.3) mas yr^ 1 . We find a relative parallax of 89 ± 18 mas in this fit, which 
is consistent with their reported value. Finally, we have determined a preliminary parallax of 
130 ± 38 mas from our own ongoing astrometry program at the Canada-France-Hawaii Telescope 
(Dupuy & Liu 2012), which further supports the unusually low luminosity of WISE J0146±4234AB. 
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3.5. Estimated Physical Properties 

We used the cloud-free Cond evolutionary models of Baraffe et al. (2003) to estimate the 
physical properties of the binary components from their luminosities and simple assumptions for 
the system age. We combine our bolometric magnitudes from above and the parallactic distance 
from Beichman et al. (2014) to compute log(Lb 0 i/^o) = —6.95 zt 0.20 dex for the primary and 
—7.01 zb 0.22 dex for the secondary. We interpolate the Cond model grid at ages of 1 Gyr and 
10 Gyr to derive masses, temperatures, and other properties from these Lboi values and present the 
results in Table 5. 

Some properties that we derive from evolutionary models vary significantly with the assumed 
age, while others do not. For example, the primary component’s mass and effective temperature 
are 4.6^}'5 Mj up and 320^4 q K at 1 Gyr but I 6 . 9 I 4 Q Mjup and 345zb45K at 10Gyr. This is because 
models predict a radius contraction of only 17% from 1 Gyr to 10 Gyr, so for a given L^oi the T e ff 
estimate changes by a small amount. Unlike the radius however, luminosity is a strong function of 
mass and age at field ages (Tbol oc M 2 4 t 1 ' 25 ; Burrows et al. 2001). Consequently, a property like 
surface gravity that depends on mass changes significantly with the assumed age (~0.7dex) while 
the effective temperature only changes by <25 K. 


3.5.1. High Surface Gravity 

There is an extensive literature on the impact of both surface gravity and metallicity variations 
on the emergent spectra of T dwarfs with T e g- = 600-1000 K (e.g., Burrows et al. 2002; Liu et al. 
2007; Leggett et al. 2007; Mace et al. 2013b), but similar studies are only beginning for <500K 
objects at the T/Y transition. The latest model atmospheres from Morley et al. (2014) show how 
the shapes of normalized near-IR flux peaks change with varying surface gravity (their Figure 15) at 
the temperatures of Y dwarfs. Models at log(g) = 5.0 dex (g in cgs units) have distinctive features 
compared to log(g) = 4.0 dex at a common T c g- = 450 K such as substantially increased flux on 
the blue side of the peak-normalized Y band, a broader J-band peak caused by increased flux on 
the red side, and a wavelength for the U-band flux peak slightly shifted to the red. Note that the 
Y-band excess is only a change in the shape of the peak and not its amplitude, as Morley et al. 
(2014) predict that Y — J actually becomes «0.5mag redder due to the integrated flux at Y band 
being more suppressed than at J band when going from log(g) = 4.0 dex to 5.0 dex. 

We observe all of these features in our integrated-light spectrum of WISE J0146+4234AB 
(Figure 2). To our knowledge, this is the first object known to possess all of these traits. A number 
of objects have been reported to display enhanced flux on the blue side of Y band (e.g., Kirkpatrick 
et al. 2011; Leggett et al. 2014a), but of these only WISE J030449.03—270508.3 (YOp) has been 
reported to show the slightly increased flux on the red side of the J-band peak (Pinfield et al. 
2014). WISE J1405+5534 (YOp) displays the shifted peak wavelength in H band (Cushing et al. 
2011) but has a normal J band and no Y-band data. Unfortunately, non-solar metallicity models 
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are not yet available at these temperatures, so it is not known if any of these traits could also be 
reproduced by a sub-solar metallicity, as is often the case for earlier type T dwarfs where the effects 
of low metallicity and high gravity can be similar (e.g., Burgasser et al. 2006a; Mace et al. 2013b). 

We conclude therefore that WISE J0146+4234AB is very likely to be relatively old, having 
higher surface gravity and/or lower metallicity compared to other known T/Y transition objects. 
The chief caveat is that we only have an integrated-light spectrum that is likely dominated by the 
primary along with our resolved broadband photometry to characterize the photospheres of the 
individual components. Our spectral decomposition analysis in Section 3.2 demonstrates that even 
a ~1 mag fainter secondary can influence the best-fit template match, since in this case it was the 
YOp WISE J1405+5534. Among the next tier of best matching templates both components of 
WISE J1217+1626AB (T8.5+Y0) were often favored, and these are also high surface gravity/low 
metallicity candidates (Liu et al. 2012; Leggett et al. 2014a). 


3.5.2. Semimajor Axis and Orbital Period 

The projected separation of WISE J0146+4234AB at discovery was p = 0 / . / 0875 ± 0 , . , 0021 
(0.93±°;i| AU), and we can use Table 6 of Dupuy & Liu (2011) to convert this into a statistical 
estimate of the semimajor axis (o). The conversion factor for very-low mass visual binaries (a/p = 
1.16jlo'3i) gives a semimajor axis estimate of (//!()/[/// (1.1 AU). However, because we discovered 
WISE J0146+4234AB very near the resolution limit of our Keck LGS AO images, we are likely 
in the “moderate discovery bias” case where the inner working angle of the discovery observations 
(IWA ~ 07055) is roughly half the size of the semimajor axis and a different conversion factor 
applies. Thus, our best estimate for the semimajor axis is 0 / . , 0951 q'q 3 ^ (I.OI 04 AU). We note that 
WISE J0146+4234AB likely has the smallest projected separation in AU at discovery of any visual 
binary ever found among L, T, or Y dwarfs. The next tightest is 2MASS J15344984—2952274AB 
that was discovered at 1.01 ± 0.03 AU (Burgasser et al. 2003). 

We can further estimate the orbital period of WISE J0146+4234AB via Kepler’s Third Law 
using the total system mass estimate we derived from evolutionary models. Assuming an age of 
10 Gyr gives M tot = 32j/g Mj up and using the smaller of the semimajor axes above (moderate 
discovery bias) yields a period estimate of P = 5.9lt§’5 yr. Even if WISE J0146+4234AB were 
much younger and thereby lower mass (M tot = 8.7l}'gMj up for an age of lGyr), its estimated 
orbital period would still be quite short (ll/jgyr). Thus, regardless of the current uncertainties in 
the system mass and semimajor axis, WISE J0146+4234AB appears to be an excellent candidate 
for dynamical mass determination via monitoring of the binary orbit over the next several years, 
as it typically requires only ~30% coverage of the orbital period to enable a direct measurement 
of the system mass (e.g., Dupuy et al. 2009a). In fact, we already detect a small amount of orbital 
motion in our Keck LGS AO imaging that spans 1.1 yr (7 ± Smasyr” 1 ). 
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4. Discussion 

Both components of WISE J0146+4234AB appear to be less luminous than any other known 
members of substellar binaries, even the Y dwarf companions CFBDSIR J1458+1013B and WISE J1217+1626B 
(Liu et al. 2011, 2012). This in turn implies that many of the properties of WISE J0146+4234AB 
are extreme relative to other known systems. To reach such low luminosities, both components of 
WISE J0146+4234AB must be planetary mass ()$13Mj up ) and/or very old. We identify several 
spectrophotometric features that are indicative of high surface gravity, so it appears likely that 
the system is indeed old, although sub-solar metallicity models at T e g = 300-500 K are needed 
to strengthen this conclusion. We estimate component masses that are quite small, with Cond 
models giving 13-21 Mj up for the primary (la range) and 12-19 Mj up for the secondary even if 
the system is 10 Gyr old. Currently, a large source of uncertainty in the component properties is 
how the integrated-light bolometric flux should be divided between the two components. Resolved 
mid-infrared photometry and near-infrared spectroscopy are needed to more precisely determine 
the physical properties of both components. 

The faintness of WISE J0146+4234A is unexpected given its T9 spectral type. Its absolute 
magnitudes are ~1.0-1.5 mag fainter in the mid-infrared and ~2mag fainter in the near-infrared 
than is typical for its spectral type, even allowing for a ±0.5 subtype uncertainty in classification. 

The faint absolute magnitudes lead to a low bolometric flux estimate, and in turn a low temperature 
of T e g = 345 ± 45 K for the primary (assuming an age of 10 Gyr). This is sslOOK colder than 
any other known late-T dwarf. 1 According to the Stefan-Boltzmann Law, at a given luminosity 
T e g- oc R 1 / 2 . Thus, a «20% discrepancy in T e g could be compensated by the radius being 40% 
smaller, i.e., 0.55 Rj up instead of 0.91 Rj up . Such a radius is unphysical as it is significantly smaller 
than any predicted or measured substellar radii (e.g., Chabrier et al. 2009; Montet et al. 2014, and 
references therein). According to evolutionary models from Saumon & Mar ley (2008), assuming 
an age older than 10 Gyr would have only a small effect ()$1%) on the radius, and assuming a 
2x larger mass would reduce the radius by only «10%. However, these effects might partially 
explain the apparent temperature discrepancy between WISE J0146±4234A and other T9/T9.5 
dwarfs. A systematic error in the Beichman et al. (2014) parallax distance could also mitigate 
this discrepancy, but to explain it entirely would require an unlikely error of a factor of >2 in the 
measured parallax. Therefore, we conclude that the discrepancy between WISE J0146±4234A’s 
spectral type and its temperature cannot be completely explained by systematic errors in the 
assumed radius or luminosity. 

One intriguing explanation for the ~100 K cooler T e g of WISE J0146±4234A compared to other 


x The lowest published Tboi-based T e ff for a late-T dwarf is 502 ± 10 K for UGPS J0722—0540 (T9), assuming an 
age of 1 Gyr (Dupuy & Kraus 2013). Examining other objects with recently published parallaxes, we compute that 
only WISE J2325—4105 (T9p), WISE J2134—7137 (T9p), and WISE J0148—7202 (T9.5) have cooler T e g estimates, 
ranging from 461 ± 13 K to 466 ± 13 K, assuming ages of 1 Gyr and using parallaxes from Tinney et al. (2014) and 
photometry from Kirkpatrick et al. (2011). 
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objects of similar spectral type is that the features in the J and H bands that are widely used to 
determine spectral types for late-T and Y dwarfs are highly dependent on gravity or metallicity. 
Since these flux peaks are shaped by molecular features (e.g., CH 4 , H 2 O, H 2 , and NH 3 ), a sensitivity 
to gravity or metallicity is not unexpected. This did not appear to be the case in the initial sample 
of late-T and Y dwarfs with Lboi-based temperatures from Dupuy & Kraus (2013) where spectral 
types tended to track well with T e g, with the possible exception of T9.5 dwarfs that seemed to be 
warmer or similar temperature compared to T9 dwarfs. However, as the sample of parallaxes for 
late-T and Y dwarfs grows, a jumbling of objects with the same T e g- but very different spectral 
types should become clearly apparent if the range of surface gravities and/or metallicities of the 
field population significantly affect near-infrared spectral types. 

With only a few of the coldest known brown dwarfs showing signatures of high surface gravity, 
it seems that objects like the components of WISE J0146+4234AB must either be relatively un¬ 
common or observationally selected against. Morley et al. (2014) predict that high surface gravity 
objects at a given temperature are indeed « 0 . 2 - 0 .3 mag fainter in the W 2 band that was used 
to discover most known objects of this type. However, this would lead to lower gravities being 
over represented by only «30%-50% in a magnitude-limited sample. In fact, given the generic 
assumption of a declining star formation history in the solar neighborhood (e.g., Aumer & Binney 
2009), older, higher gravity brown dwarfs should be much more common at a given T e g than lower 
gravity ones, unless the mass function compensates by having a strong preference for producing 
lower mass brown dwarfs. If higher gravity brown dwarfs are uncommon it would imply either an 
unusual age distribution for brown dwarfs in the solar neighborhood (favoring younger ages) or a 
large undiscovered population of planetary-mass brown dwarfs that are the older and colder but 
similar mass counterparts to the current Y dwarf sample. 

Fortunately, WISE J0146+4234AB will soon provide an opportunity to much better constrain 
its physical properties. With a projected separation of 0.93]/g'}g AU and estimated orbital period 
<10 yr, WISE JO 146+4234AB is not only a contender for the tightest known ultracool binary, it 
is also the first binary discovery at the T/Y transition that is likely to yield a dynamical mass 
quickly, as other discoveries have estimated orbital periods of ~40-400yr. The relative orbit of 
WISE J0146+4234AB could be determined by astrometric monitoring within just a few years. A 
dynamical mass for WISE J0146+4234AB will test the hypothesis that the components are old and 
massive with high surface gravity, which is solely based on current theoretical predictions of the 
spectrophotometric properties of brown dwarfs. As perhaps the only planetary-mass brown dwarfs 
with a dynamical mass measurement for the next several years, WISE J0146+4234AB will be a 
prime target for James Webb Space Telescope spectroscopy to test cold model atmospheres where 
the physical parameters are well constrained by observations. 
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Fig. 1.— Contour plots of our Keck LGS AO images from which we derive astrometry and flux 
ratios (Table 1). Contours are in logarithmic intervals from unity to 10% of the peak flux in each 
band. Images have been rotated such that north is up. 
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Fig. 2.— Top: Integrated-light spectrum of WISE J0146+4234AB obtained with Gemini/GNIRS 
(smoothed to R ~ 150; black). The spectrum is flux calibrated using our Gemini/NIRI photom¬ 
etry (Table 2). Plotted for comparison are spectral standards that have been normalized to the 
peak flux of WISE J0146+4234AB at 1.15-1.40/rm. Bottom: Zoomed in plots of our spectrum 
compared to the same standards, where all spectra are instead normalized to their peak in that 
band. Differences between our observed spectrum and the standards are shown with error bars in 
the bottom subpanels. We plot one data point per resolution element. 
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Fig. 3.— Integrated-light spectrum (black) and best matching component templates (colored lines). 
The bottom subpanel shows the observed J-, H-, and CH^s-band broadband flux ratios used to 
constrain the decomposition (filled black circles with errors) and the resulting flux ratios computed 
from the best matching template pair (open brown squares). 
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Fig. 4.— Color-magnitude diagrams for all objects with spectral types T8 and later that have ac¬ 
curate photometry (Leggett et al. 2014b, and references therein) and direct distance measurements 
(Beichman et al. 2014; Dupuy & Kraus 2013; Dupuy &; Liu 2012; Leggett et al. 2012; Luhman 
et al. 2011; Luhman & Esplin 2014; Tinney et al. 2014; Wright et al. 2013). Binary systems are 
shown in both integrated light and as resolved components when possible. WISE J0146+4234AB 
and its components are plotted as star symbols, showing unusual absolute magnitudes (and some¬ 
times colors) for their spectral types. Other binaries are plotted as diamonds. All data points 
are color coded according to spectral type, with open/white points indicating that no spectra are 
available. Near-infrared photometry is on the Mauna Kea Observatory (MKO) system, and objects 
with distance modulus errors >0.5 mag are not plotted for clarity. 
















Table 1. Keck LGS AO Observations of WISE J0146+4234AB 


Date 

(UT) 


Airmass 

Filter 

FWHM 

(mas) 

P 

(mas) 

P.A. 

(°) 

Am 

(mag) 

2012 Sep 

7 

1.084 

Y 

48 ±26 

85.9 ± 4.4 

265.0 ±4.4 

0.83 ±0.23 

2012 Sep 

7 

1.118 

J 

58 ± 14 

88.2 ±3.4 

258.8 ±2.0 

0.98 ± 0.03 

2012 Sep 

7 

1.091 

H 

52 ±6 

87.9 ±3.5 

259.5 ± 1.8 

0.97 ±0.10 

2012 Sep 

7 


mean 


87.5 ±2.1 

259.7 ± 1.3 


2012 Oct 

8 

1.093 

Y 

64 ± 17 

91.0 ±2.9 

260.5 ±4.1 

0.83 ±0.23 

2012 Oct 

8 

1.236 

J 

60 ± 16 

93.9 ±4.2 

262 ±6 

1.07 ±0.08 

2012 Oct 

8 

1.198 

H 

49 ±4 

85.6 ±3.6 

258.3 ±4.1 

1.15 ±0.14 

2012 Oct 

8 

1.129 

CH 4 s 

55 ±11 

93.5 ± 4.4 

262.0 ± 1.2 

1.32 ±0.13 

2012 Oct 

8 


mean 


90.6 ±1.8 

261.6 ±1.1 


2013 Oct 

22 

1.098 

j 

56 ± 10 

92.9 ±4.0 

263.6 ± 1.3 

0.90 ±0.16 


Note. - All photometry on the MKO photometric system. The identical K-band 
flux ratios are not a typo. 
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Table 2. Gemini-North NIRI Photometry of WISE J0146+4234 


Filter 

^exp 

(s) 

Photometry 

(mag) 

Date 

(UT) 

Y 

6 x 60 

21.60 ±0.15 

2013 Jan 10 

J 

6 x 60 

20.69 ± 0.07 

2013 Jan 10 

H 

28 x 30 

21.30 ±0.12 

2013 Jan 10 

CH 4 s 

28 x 30 

20.51 ±0.14 

2013 Jan 12 

K 

28 x 30 

22.4 ± 0.4 a 

2013 Jan 10 


a Given the very low signal-to-noise ratio at K 
band, we suspect our photometric errors in this 
filter may be underestimated. 


Table 3. Median-Flux Spectral Indices for WISE J0146+4234AB 


Spectral 

Index 

WISE J0146+4234AB 

Average Median-Flux Values from Mace et al. (2013a) 

Best-match 
Spec. Type 

T8.5 

T9 

T9.5 

Y0 

Wj 

0.188 ±0.031 

0.295 ±0.047 

0.203 ± 0.038 

0.145 ±0.021 

0.117 ±0.043 

T9 

J-narrow 

0.923 ± 0.036 

0.884 ± 0.044 

0.879 ± 0.053 

0.865 ± 0.071 

0.778 ± 0.050 

T9 

CH 4 - J 

0.063 ± 0.024 

0.188 ±0.056 

0.117 ±0.030 

0.071 ±0.032 

0.045 ± 0.031 

T9.5 

ch 4 -h 

0.120 ±0.029 

0.122 ±0.035 

0.068 ± 0.067 

-0.021 ±0.116 

0.066 ±0.133 

T8.5 

NH a-H 

0.457 ±0.055 

0.610 ±0.059 

0.539 ± 0.054 

0.443 ±0.101 

0.385 ±0.105 

T9.5 

Y/J 

0.571 ±0.038 

0.432 ± 0.082 

0.448 ± 0.095 

0.357 ±0.008 

0.423 ±0.131 


H/J 

0.539 ±0.027 

0.503 ± 0.045 

0.555 ±0.030 

0.504 ± 0.080 

0.467 ± 0.061 


K/J 

0.066 ± 0.007 

0.116 ±0.039 

0.100 ±0.032 

0.033 ± 0.083 

0.063 ± 0.031 


J- wing 

0.322 ± 0.043 


0.320 ± 0.005 a 


0.164 ±0.012 a 

T9 


Note. — All spectral indices except J-wing are described in Mace et al. (2013a). We exclude the indices H 2 O — J 
and H 2 O — H as they are both saturated, i.e., have zero flux in the numerator, for WISE J0146+4234AB. 

a These are the values for the spectral standards UGPS J0722—0540 (T9) and WISE J1738+2732 (Y0) reported 
by Pinfield et al. (2014) who originally defined this index. 
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Table 4. Integrated-Light Properties of WISE J0146+4234AB 


Property 

A+B 

Ref. 

Measured 

33DSS (mag) 

24.10 ±0.13 

Lodl3 

Lmko (mag) 

21.60 ±0.15 

§2.2 

Jmko (mag) 

20.69 ± 0.07 

§2.2 

Hmko (mag) 

21.30 ±0.12 

§2.2 

CH 4 smko (mag) 

20.51 ±0.14 

§2.2 

A'mko (mag) 

21.75 ±0.25 

§2.3 

[3.6]irac (mag) 

17.50 ±0.07 

Kir 12 

[4.5]irac (mag) 

15.08 ±0.02 

Kir 12 

Spectral type 

T9p 

§2.3, §3.2 

Distance (pc) 

10-6111 

Beil4 

IP* cos 5 ("yr” 1 ) 

-0.441 ± 0.013 

Beil4 

vs ("yr 1 ) 

-0.026 ±0.016 

Beil4 


Estimated 

m h oi (mag) 
log(±bol/ Lq) 

21.58 ±0.12 

—6 fi7+ 0 - 12 

D - D ‘ -0.15 

§3.4 

§3.4, Beil4 


References. — References: § numbers refer to this 
paper; Beil4 (Beichman et al. 2014); Kir 12 (Kirk¬ 
patrick et al. 2012); Lodl3 (Lodieu et al. 2013). 
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Table 5. Resolved Properties of WISE J0146+4234AB 


Property 

Component A 


Component B 

Measured 

Separation * 1 

0"0875 ± 

0"0021 

P.A. a 

259?7 ± 

1?3 

AImko (mag) 

0.83 ± 

0.16 b 

AJmko (mag) 

0.99 ± 

0.03 b 

ARmko (mag) 

1.03 ± 

0.08 b 

A.CH 4 SMKO (mag) 

1.32 ± 

0.13 

Pmko (mag) 

22.02 ±0.16 


22.85 ±0.19 

■Imko (mag) 

21.06 ±0.07 


22.05 ±0.07 

Rmko (mag) 

21.66 ± 0.12 


22.69 ±0.14 

CH 4 SMK 0 (mag) 

20.79 ±0.14 


22.11 ±0.17 

Spectral type 

T9 


Y0 

Estimated 

Am bo i (mag) 

0.2 ± 0.8 

m bo i (mag) 

22.2 ±0.4 


22.4 ±0.4 

log(L bo i/L 0 ) 

-6.95 ±0.20 


-7.01 ±0.22 

Model-derived (Cond, t = 

: 1 Gyr) 

Mtot (Mj up ) 

0 '7+1.3 

'- 1.6 

Mass (Mj up ) 

4.6±i;° 


4.3ll° 

T eS (K) 

3201^5 


3 IOI 40 

Radius (i?j up ) 

1 f)674 +0 ' 0027 
i.uo /+_ 0 0017 


1.0675lg;ooil 

log (g) (cgs) 

3 QQ + 0 - 11 
°‘ yy 0.09 


3-961°+ 

Model-derived (Cond, t = 

10 Gyr) 

M to t (Afj U p) 

32i^ 

Mass (M Jup ) 

16.911® 


15.9±|;t 

T eS (K) 

345 ± 45 


330 ± 45 

Radius (i?j U p) 

0.913±8;811 


f) Q 1 0+0.034 

u.yiy_o.oi5 

log(s) (cgs) 

4.69lo.ii 


4 6 c ) +0 14 
^• oo -0.12 


a Weighted average of measurements in different filters 
at epoch 2012 Sep 7 UT. 

b Weighted average of measurements from multiple 
epochs. 













